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Irradiation with Pb ions was used to study the efTect of disorder on the in-plane London penetra- 
tion depth, A(r), in single crystals of Ba(Fei_a;Ta;)2As2 (T=Co, Ni). An increase of the irradiation 
dose results in a monotonic decrease of the superconducting transition temperature, T^, without af- 
fecting much the transition width. In both Co and Ni doped systems we find a power-law behavior, 
AA(r) oc T", with the exponent n systematically decreasing with the increase of disorder. This 
observation, supported by the theoretical analysis, conclusively points to a nodeless state with 
pairbreaking impurity scattering (interband) with strength being intermediate between Born and 
unitary limits. 



The mechanism of superconductivity in Fe-based su- 
perconductors |l| is a focus of extensive experimental 
Q and theoretical 0, |^ efforts. Understanding the su- 
perconducting gap structure is crucial for identifying the 
pairing mechanism. While most of the experiments sug- 
gest that the general structure of the pairing state be- 
longs to the most symmetric Aig class, the multi-band 
nature of these materials allows for a number of possibil- 
ities, including the conventional s-wave state, extended 
s^-state with different signs of the order parameter on 
two bands 0-0] and states with highly anisotropic gaps 
and even nodes [S- 12|. 

The London penetration depth can be measured with 
great precision and its variation with temperature de- 
pends sensitively on the gap structure. For T < Tc/3, 
a conventional isotropic s-wave gap Aq results in an ex- 
ponential behavior, AA(r) cx exp(— Aq/T), which is pre- 
served even with the addition of non-magnetic impurities 
[isj . Unconventional pairing states, on the other hand, 
are susceptible to the presence of non-magnetic impu- 
rities [3, lj-[l6|. In nodal d-wave superconductors, A(r) 
exhibits a power-law behavior, AA(T) c>c T", with the ex- 
ponent n increasing from n=\ in the clean case to n—2 
in the dirty limit [17|. For the extended state, the 
opposite trend is expected: AA(T) is exponential in the 
clean limit, changing with disorder to a power-law, with 
n as low as 1.6 |14l. Il5|. 

Experimentally, a power-law behavior with the expo- 
nent 2 < n < 3 has been observed in most of the iron- 
based superconductors [lS'-^25!j. This characteristic ex- 
ponent n ^2 can be explained in both dirty d-wave and 
dirty scenarios. The disorder is always present in 
the iron-based superconductors where doping (e.g., Co or 
Ni in this work) is needed to induce superconductivity. 
One way to resolve this complication is to deliberately in- 
troduce defects that do not contribute extra charge but 
rather only increase the scattering rate. Various ways 
of controlling the scattering rate have been suggested in 



earlier studies, especially for the high-Tc cuprates, and 
the effects have been examined by using transport and 
magnetic measurements [2^, 113. In particular, irradia- 
tion with heavy ions, besides producing efficient pinning 
centers, also significantly enhances scattering, as evident 
from the significant increase of the normal state resistiv- 
ity H, 2^ as well as the suppression of Tc. 

In this work, we study the in-plane London penetra- 
tion depth in single crystals of optimally Co-and Ni- 
doped BaFe2As2 superconductors irradiated with 1.4 
GeV 208p|-,56-i- ions. In both systems we find mono- 
tonic suppression of Tc with the increase of the irradi- 
ation dose without notable transition broadening. The 
London penetration depth exhibits a power-law behav- 
ior, AA(T) cx AT" (2.2 < n < 2.8), with the exponent 
n decreasing with the irradiation dose. This observa- 
tion, supported by the theoretical analysis, provides the 
most convincing case for the nodeless state with pair- 
breaking scattering (interband) of intermediate strength, 
between Born and unitary scattering limits. 

Single crystals of Ba(Fei_^T^)2As2 (T=Co, Ni de- 
noted FeCol22 and FeNil22, respectively) were grown 
out of FeAs flux using a high temperature solution growth 
technique 3Q , i31i] . X-ray diffraction, resistivity, magneti- 
zation and wavelength dispersive spectroscopy (WDS) el- 
emental analysis have all shown good quality single crys- 
tals at the optimal dopings with a small variation of the 
dopant concentration over the sample and sharp super- 
conducting transitions, Tc= 22.5 K for FeCol22 and 18.9 
K for FeNil22 [s^, 31 1. The in-plane London penetra- 



tion depth was measured by using the tunnel diode res- 
onator technique [3^-[35|. The sample was placed with 
its crystallographic c-axis parallel to a small excitation 
field, Hac ~ 20 mOe. The shift of the resonant frequency, 
Af(T), is proportional to the sample magnetic suscepti- 
bility, xiT) via A/(T) = -GAttxIt). Here G is a geo- 
metric calibration factor, G — foVs/2Vc{l~ N), where N 
is the demagnetization factor, Vs is the sample volume, 



2 



and Vc is the coil volume. The calibration factor was 
determined from the full frequency change by physically 
pulling the sample out of the coil. The magnetic suscep- 
tibility can be written in terms of A and the characteristic 
length R, Attx =_(A/i?) tanh(i?/A) - 1, from which AA 
can be acquired [33| . 

To examine the effect of irradiation, ~ 2 x 0.5 x 0.02 — 
0.05 mm^ single crystals were selected and then cut into 
several pieces preserving the width and the thickness. We 
compare sets of samples, where the samples in each set 
are parts of the same original large crystal. Several such 
sets were prepared and a reference piece was kept unirra- 
diated from each set. Irradiation with 1.4 GeV 208p|356+ 
ions was performed at the Argonne Tandem Linear Ac- 
celerator System (ATLAS) with an ion flux of ~ 5 x 10^^ 
ions-s~"'^-m~^. The actual total dose was recorded in each 
run. The density of defects (c?) created by the irradia- 
tion is usually expressed in terms of the matching field, 
Btf, = ^od, which is obtained assuming one flux quanta, 
$0 ~ 2.07 X 10~^ G-cm^ per ion track. Here we studied 
samples with = 0.5, 1.0 and 2.0 T corresponding to 



^. The sample thickness was chosen in the range 
20 — 50/xm to be smaller than the ion penetration 



cm 
of 

depth, 60 — 70 /im. The same samples were studied by 
magneto-optical imaging. The strong Meissner screening 
and large uniform enhancement of pinning have shown 
that the irradiation has produced uniformly distributed 
defects [36|. 

Figure [T] shows AA(T) for FeCol22 (top panel) and 
FeNil22 (bottom panel). The low-temperature region 
up to ss Tc/3 is shown in the main frame of each panel. 
Vertical offsets were applied for clarity. The normalized 
penetration depths in the vicinity of Tc are shown in the 
inset of each panel to highlight the suppression of Tc as 
the radiation dose increases. Whereas Tc is clearly sup- 
pressed, the transition width remains nearly the same 
(see Fig. [3] below). All samples exhibit a power-law vari- 
ation of AA(T) cx T" with 2.5 < n < 2.8 up to T^/S, 
while the exponential fitting failed in all cases. The 
best fitting curves are shown by solid lines in Fig. [1] 
We note that the present set of FeCol22 samples ex- 
hibits higher exponents, n, compared to previous works 
This variation of n is likely due to disorder, 
as we clearly demonstrate in this work. Consequently, 
it is important to conduct a comparison of radiation ef- 
fects on the samples made from the same large crystal. 
Magneto-optical characterization has shown a homoge- 
neous superconducting response [s^ and the widths of 
the superconducting transitions were much smaller than 
the absolute shift due to irradiation, see Fig. |31 There- 
fore, it is very likely that the effects reported here are 
caused by the enhanced scattering induced by the heavy 
ion bombardment. 

To further analyze the power-law behavior and its 
variation with irradiation, we plot AA as a function of 
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FIG. 1. Variation of the in-plane London penetration depth, 
AA(r), for irradiated FeCol22 (top panel) and FeNil22 (bot- 
tom panel). The low-temperature variations are shown in the 
main frame of each panel along with the best power-law fits. 
The curves are offset vertically for clarity. The variations in 
the vicinity of Tc are shown in the insets of each panel. 



{t = T/Tcy-" in Fig.H where the riQ values for FeCol22 
and FeNiI22 were chosen from the best power-law fits of 
the unirradiated samples (see Fig.[3l). While the data for 
unirradiated samples appear as almost perfect straight 
lines showing robust power-law behavior, the curves for 
irradiated samples show downturns at low temperatures 
indicating smaller exponents. This observation, empha- 
sized by the plots of the derivatives dAX{t)/dt"''' in the 
inset of Fig. [21 points to a significant change in the low- 
energy excitations with radiation. 

The variations of Tc and n upon irradiation are illus- 
trated in Fig. [21 Dashed lines and circles show FeCoI22, 
while solid lines and triangles show FeNiI22. The up- 
per panel shows the variation of Tc and the width of 
the transition. Since i?^ is directly proportional to the 
area density of the ions, d, we can say that Tc decreases 
roughly linearly with d. The same trend is evident for 
the exponent n shown in the lower panel of Fig. [3l The 
fitting pre-factor A increases somewhat upon the increase 
of irradiation dose, but remains smaller than the value 
measured previously in unirradiated samples fisl [Tol [13] . 

The experimental results fit comfortably within the hy- 
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FIG. 2. Detailed comparison of the functional form of AA(r) 
for irradiated FeCol22 and FeNil22. In the main panels 
AA(r) is plotted vs. {t = T/Tc)"° with the exponents n taken 
from the best fits of unirradiated samples: no = 2.8 and 2.5 
for FeCol22 and FeNil22, respectively (see Fig[l]). Appar- 
ently, irradiation causes low-temperature deviations, which 
are better seen in the derivatives, dAX{t)/dt"° , plotted in the 
insets. 



pothesis of superconductivity with two isotropic gaps. 
The superfluid density in linear response is 



'N 



i=l,2 



#[v/,j(g)V/,^]a 



■(e^+Af)3/2 



(1) 

where one sums over the contributions from the elec- 
tron and hole bands; vj^^ and Nf^i are the Fermi ve- 
locity and density of states in these bands, taken to 
be equal for the calculations. Two order parameters 
Ai,2 are computed self-consistently together with the t- 
matrix treatment of impurity effects, which renormalize 
the Matsubara energies iim = iSm — '^imp,i and the gaps 
A,; = A,;-f Aimp,i[37|. Impurities are characterized by the 



strength of the potential for scattering within each band, 
vii{= V22), given by the phase shift 5 — tan~^(7rA^/-yii), 
the ratio of potentials for inter-band and intra-band hop- 
ping, 5v = wi2/i'ii, and the impurity scattering rate 




FIG. 3. Top panel: The suppression of Tc with disorder rela- 
tive to unirradiated T"". The vertical bars denote the width 
of the transition and correspond to temperatures where the 
diamagnetic response changed from 90% (onset) to 20% (end 
of the transition), see insets in Fig.[T] Lower panel: exponent 
n vs. B^. 



The essential theoretical results are presented in Fig.|3] 
while the full calculations will be published elsewhere 
(ssj . The best agreement with the experiment is obtained 
for two isotropic gaps, A2 ~ — O.6A1, strong inter-band 
scattering 6v — 0.9 and phase shift S = 60° between 
the Born {S -> 0) and unitary limits {6 -> 90°). The 
calculated p{T) was fitted to the power-law, p{T)/po « 
p{0)/po~a (T/Tco)", which is directly related to the pen- 
etration depth, AA(T)/Ao « a'(T/Tc)", with po and Aq 
being the T = superfluid density and penetration depth 
in the clean system and a' = {a/2)[Tc/Tco]"[po/ piO)]^^^ ■ 
We find that with an increase in F, the power n decreases 
from n > 3 to n « 2 (see Fig. H^a)), which is in per- 
fect agreement with experiment. The values of n depend 
sensitively on the structure of the low-energy density of 
states, which is shown in Fig. |H[b). The intermediate 
strength of scatterers is important for creation of a small 
band of mid-gap states separated from the continuum. 
As the disorder is increased, these states close the gap 
in the spectrum and gradually increase in magnitude, 
driving the low-temperature power-law dependence from 
exponential-Hke, n > 3, to n w 2. This behavior of n is 
largely independent of the details of the model, whereas 
a' can slightly increase or decrease depending on the dif- 
ferent ratios of the gaps on two Fermi surfaces and dif- 
ferent impurity parameters. 

Finally, in Fig. ^c) we show the central result of our 
study: the correlation between Tc and n. Note that these 
two quantities are obtained essentially independently of 
each other. Assuming that the unirradiated samples have 
some disorder due to doping, and scaling T"" to lie on 
the theoretical curve, we find that the Tc{Bij,) of the irra- 
diated samples also follows this curve. The assumption 
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FIG. 4. (a) Superfluid density and (b) the density of states 
n{e) — N{e)/Nf, computed for the state with sign- 
changing isotropic gaps and strong interband impurity scat- 
tering, between the Born and unitary hmits. The dashed 
line in (a) is an example of a power-law fit p{T)/po — 
p(0)/po - a{T/Tco)" for < r < 0.4 Tco; best fitting param- 
eters for a given set of F = Uimp/irNf are listed in the table, 
(b) As the impurity concentration riimp ~ T)) increases, the 
band of mid-gap states approaches the Fermi level and the ex- 
ponent n is reduced, (c) Tc vs. power n, from the theoretical 
model (triangles) and experiment (squares and circles). 



of similarity between doping and radiation-induced disor- 
der, implied in this comparison, while not unreasonable, 
deserves further scrutiny. 

In summary, we determined the effect of irradiation on 
A(T) and demonstrated that the disorder-induced reduc- 
tion of the power-law exponent 2 < n < 3 is naturally 
explained in terms of the isotropic extended s-wave state 
[3, 01 with pair- breaking interband scattering!^ IigI. Isoj. 
We have also considered models for nodal states, but they 
showed the opposite trends: increase of n with disorder in 
the interval 1 < n < 2, and thus can be excluded. Taken 
together with reports of fully gapped states from ther- 
mal conductivity j40i] and angle resolved photoemission 
spectroscopy [4l|, our results present the most convinc- 



ing case, up to date, in favor of the extended s pairing 
symmetry having a nodeless order parameter in the op- 
timally doped 122 system. 

The picture of strong pair-breaking scattering is also 
consistent with recent proposals of the universal be- 
havior in the thermal and electromagnetic responses of 
iron-based superconductors [H, m, Nonetheless, 
we should note that nodal states may still exist in P- 
containing compounds [44 46] or along the c-axis of heav- 
ily overdoped Bal22 pnictides [23j . 
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